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ABSTRACT 

Transiting planets, such as HD 209458b, offer a unique opportunity to scrutinize the planetary atmospheric content. 
Aithough molecular hydrogen is expected to be the main atmospheric constituent, H2 remains uncovered because of the 
lack of strong transition from near-ultraviolet to near-infrared. Here we analyse the absorption spectrum of HD 209458b 
obtained by Sing et al. (2008a) which provides a measurement of the absorption depth in the 3000-6200 A wavelength 
range. We show that the rise in absorption depth at short wavelengths can be interpreted as Rayleigh scattering 
within the atmosphere of HD 209458b. Since Rayleigh scattering traces the entire atmosphere, this detection enables 
a direct determination of the pressure-altitude relationship, which is required to determine the absolute fraction of 
other elements such as sodium. At the zero altitude defined by the absorption depth of 1.453%, which corresponds 
to a planetary radius of 0.1205 times the stellar radius, we find a pressure of 33±5mbar. Using the variation of the 
Rayleigh scattering cross-section as a function of wavelength, we determine the temperature to be 2200±260 K at 
33mbar pressure. 
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1. Introduction 

Transiting planets such as HD 209458b provide a unique 
opportunity to scrutinize their atmospheric content {e.g., 
'. Charbonneau et al. 2002, Vidal-Madjar et al. 2003, 2004, 
2008; Ballester et al. 2007; Barman 2007). The theory 
\ of transmission spectroscopy was developed in detail by 
. the pioneering works of Saeger & Sasselov (2000), Brown 
' (2001), and Hubbard et al. (2001). Primary transits can re- 
I veal minute quantities of gas through spectral absorption 

■ leading to larger apparent planet size, or larger absorption 
' depth in the transit light curve, at wavelengths characteris- 
, tic of the absorber. Due to Rayleigh scattering, the main at- 

■ mospheric constituent, H2, is seen by means of the ramp-up 
I of apparent planetary radius from 500 to 300 nm (Hubbard 

et al. 2001). Using transmission spectroscopy, several mi- 
nor constituents of the atmosphere of HD 209458b have al- 
ready been uncovered, including sodium, atomic hydrogen, 
carbon and oxygen. 

The extrasolar planet HD 209458b was observed during 
primary transit using the Hubble Space Telescope STIS 
instrument, at low resolving power using the G430L and 
G750L spectrograph, and at medium resolution using the 
G750M spectrograph. The data was analyzed as described 
in Balleter et al. (2006) and Sing et al. (2008a); special care 
was taken to correct for the limb-darkening effect. In this ar- 
ticle, we use the absorption depth (AD) of the planet, which 
was measured as a function of wavelength from 3000 A to 
-6200 A. 

In the whole paper, the zero altitude in the planet at- 
mosphere is defined to be the mean absorption depth of 
1.453%, measured by Knutson et al. (2007). This corre- 
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sponds to a planetary radius of 0.1205 times the stellar 
radius. 

In Sect [21 we present first estimates of temperature and 
pressure that can be obtained assuming that the Rayleigh 
scattering by molecular hydrogen is responsible for the ob- 
served increase in absorption depth at short wavelengths. 
In Sect 121 we present a fit to the data including Nal ab- 
sorption lines. Discussion and conclusion can be found in 
Sect. 1 and 

2. First estimates 

2.1. The HD 209458b transit spectrum from 3000 to 6200 A 

In a plot of the absorption depth as a function of wave- 
length, an increase in absorption depth is observed to- 
ward shortest wavelengths below '^5000 A (Fig. IT]). Above 
5000 A, the absorption depth remains roughly constant up 
to 5500 A, before rising again from 5500 to 6000 A. This 
last increase in absorption depth, which has a maximum at 
about 5900 A is interpreted to be due to a high abundance 
of sodium in the bottom atmosphere, below the Na2S cloud 
level (Sing et al. 2008b). 

In the short wavelength side, the mean absorption 
depth is measured to be 1.444%-h/-0.002% between 4500 
and 5500 A, while it is 1.475%±0.005% between 3000 and 
3900 A. This significant variation of 0.030% in absorption 
depth was already identified in the same dataset (Ballester 
et al. 2007). Ballester et al. (2007) proposed that the in- 
crease in absorption depth below 4OOO A could be due to 
absorption in the Balmer series by an optically thin layer 
of excited hydrogen atoms at high altitude in the plan- 
etary atmosphere. Here we propose another explanation: 
Rayleigh scattering. 
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2.2. Rayleigh scattering 

Following the derivation of Lecavelier des Etangs et al. 
(2008), we define Tcq as the optical depth at the altitude 
Zeq, such that a sharp occulting disk of radius Rp + Zcq pro- 
duces the same absorption depth (AD) as the planet with a 
translucent atmosphere. In other words, Teq is defined to be 
AD = {Rp + z{T = req))^/i?stari whcrc Rp and i?star are the 
planet and star radii, respectively. Numerical integration of 
absorption depth from model atmospheres shows that Tcq is 
quasi-constant; its value is ~0.56 for a wide range of atmo- 
spheric scale height, H, provided Rp/H is between ~30 and 
~3000 (Lecavelier des Etangs et al. 2008). For a tempera- 
ture T, the atmosphere scale height is given hy H = kT/ fj,g, 
where g is the gravity and /i is the mean molecular mass. 
In the case of HD 209458b, Rp/H varies between 100 and 
1000, when the temperature varies from 300 to 3000 K. For 
a given atmospheric structure and composition, the theo- 
retical absorption depth at a wavelength A can be therefore 
calculated by finding z{t = Tcq, A), which solves the equa- 
tion t(z. A) = Tcq = 0.56 (for details, see Lecavelier et al. 
2008). 

The Rayleigh scattering cross-section follows a power- 
law function of wavelength to the power of four, of the 
form: (t(A) = ao{X/Xo)~'^, where ctq is the Rayleigh scat- 
tering cross-section at a reference wavelength Aq. The high 
value of exponent in the power law is the reason behind the 
Earth's blue sky. In this case, equations from Lecavelier des 
Etangs et al. (2008) that give the planet radius as a func- 
tion of the wavelength can be solved analytically. For the 
absorption depth of Rayleigh scattering, this leads to the 
simple equation 



AD = ADoi I -^In^ 



Rr, 



Ao 



(1) 



where ADq is the absorption depth at a reference wave- 
length Aq. 

Two important quantities appear in this equation. First, 
the absorption depth follows a linear relation as a function 
of the logarithm of wavelength; the slope is characteris- 
tic of the scale height H, which is directly proportional to 
the temperature. Second, the total density and therefore 
the pressure at the altitude defined by ADq corresponds to 
density and pressure at which Rayleigh scattering is opaque 
at wavelength Aq. Thus, the determination of the couple 
{ADq,Xq) provides an absolute reference for the altitude- 
pressure relationship. This absolute reference is required 
to interpret measurements of relative absorption depth by 
minor species such as sodium, and to determine their abun- 
dances (Sing et al. 2008b). 

The temperature is directly related to the observed vari- 
ation in absorption depth as a function of wavelength. For 
HD 209458b, assuming a temperature of 1500 K that pro- 
vides a typical scale height of 550 km with the planetary 
parameters of Knutson et al. (2007), and using a mean 
molecular mass of fi—2.3'mp, where rup is the mass of the 
proton, the variation in absorption depth is expected to 
be to be -0.034% from 3000 A to 5000 A; this is similar to 
what is observed in the G430L data. In the other way, if 
interpreted in term of Rayleigh scattering, the observed in- 
crease in absorption depth at short wavelengths of 0.031% 
between 3500 and 5000 A (SectlT]) implies a temperature of 
about 2000 K. 



Concerning the pressure, the fit of the measured ab- 
sorption depth from 3200 to 5000 A by the law given 
in Eq. [1] provides ADq =1.453% at Ao=4300A, where 
1.453% corresponds to the mean planet radius measured 
by Knutson et al. (2007). This allows the derivation of 
the pressure at the planet radius (zero altitude) to be 
Pq — Tcq/(JQ X y/kTfj,g/2TTRp, where ao is the Rayleigh 
scattering cross-section at 4300 A (Lecavelier des Etangs et 
al. 2008). Using the refractive index of molecular hydrogen 
(r - 1) = 1.32 X lO-"*, one derives ctq = 2.3 x lO^^^cm^ at 
Ao=4300A, and Pq ~30 millibars. 

3. Fit of the spectrum 

3.1. Fit witli Rayleigh scattering and sodium line 
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Fig. 1. Plot of absorption depth as a function of wave- 
length (histogram). The data were rebinned by 20 pix- 
els corresponding to 55 A per bin at wavelengths below 
5500 A. The measured absorption depths have typical 1- 
cr error bars of 0.015% to 0.020% per bin below 4000 A, 
and about 0.010% above 4500 A. The thick line shows the 
best fit to the data using the temperature-pressure-altitude 
profile given in Sing et al. (2008b). The dashed line shows 
the absorption with Rayleigh scattering only. 

Beyond simple estimates using a model that includes 
only Rayleigh scattering, we can measure the atmospheric 
structure of HD 209458b more accurately using a global fit 
to the data which introduces other absorbents, and includes 
medium-resolution data which further constrains the char- 
acteristics of sodium in the high atmosphere (Sing et al. 
2008b). With other absorbents, the absorption depth can- 
not be calculated analytically and a numerical fit to the 
data is required. 

We fit the whole set of data, including data acquired 
using the medium resolution G750M spectrograph (Sing 
et al. 2008a, 2008b). In the fit, the pressure and tempera- 
ture profile as a function of altitude, as well as the sodium 
abundance, are free parameters. There is a total of 7 free 
parameters that are described in Sing et al. (2008b): five pa- 
rameters describe the T-P profile and two parameters, the 
sodium abundance below and above the sodium conden- 
sation temperature. The sodium absorption line profile is 
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calculated using coUisional broadening (Iro et al. 2005). We 
obtain a satisfactory fit, which has a of 40 for 52 degrees 
of freedom in the low- resolution data set (Fig. [T]) . The re- 
sults for sodium clouds and the corresponding temperature- 
pressure diagram are given and discussed in Sing et al. 
(2008b). 

Defining zero altitude to be an absorption depth of 
1.453%, we find that the temperature and pressure at 
this altitude are constrained mainly by Rayleigh scatter- 
ing below 5000 A. We find Pq = 33 ± 5 millibars and 
To = 2200 ± 260 K (l-a). At this pressure and tempera- 
ture, the Rayleigh scattering makes the atmosphere opti- 
cally thick at wavelength shorter than '^5000 A, and causes 
a significant increase in the absorption depth of '^0.03% 
from 5000 A to 3000 A (Fig.[T|). The blue wing of the sodium 
line overcomes the absorption due to Rayleigh scattering at 
wavelengths longer than ^5000 A. An additional absorption 
in the red wing of the sodium line, beyond 6200 A, could be 
explained by TiO and VO molecules (Desert et al., in prepa- 
ration). The sodium abundance required to fit the blue tail 
of the line profile is found to be A:(NaI) = 4.3li;3 ^ 10~^ 
which is two to four times the solar abundance (Asplund el 
al. 2005). 

Error bars for temperature and pressure at the zero al- 
titude level, are plotted in FigUl These error bars are calcu- 
lated using the difference between the data and a fit to 
the absorption depth measurements from 3200 to 6200 A. 
For a given temperature and pressure at zero altitude, the 
other parameters of the fit are free to vary. Therefore, the 
resulting error bars include the uncertainties in the sodium 
abundance and other parameters of the T-P profile. 
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Fig. 2. Plot of error bars for the temperature and the pres- 
sure at the zero altitude level corresponding to the altitude 
at which the absorption depth is 1.453%. 

This fit to the data does not assume an isothermal tem- 
perature profile as function of altitude in the atmosphere. 
Even more, the temperature gradient is found to be large 
and close to the limit allowed by the adiabatic gradient 
(Sing et al. 2008b). An isothermal assumption is used only 
to calculate the column density in the limb grazing line of 
sight. Following Fortney (2005), we used an isothermal de- 
crease of density with altitude to derive an approximation of 
the column density, Niso = n^^-nHKp^ where n is the vol- 



ume density at lowest altitude in the line of sight (see also 
Lecavelier des Etangs et al. 2008). With a temperature gra- 
dient dT / dz, a better approximation of the column density 
is given by iV « N-.^o ^jTdz/HdT F ( 1 /2 Tdz /HdT) /F ( 1 + 
Tdz/HdT), that is iV « iViso(l - 0.125HdT/Tdz). For the 
maximum temperature gradient that corresponds to the on- 
set of convection (adiabatic T-P profile) and is given by 
HdT/Tdz = 2/7, the difference between A'iso and the ac- 
tual column density is less than 4%. The results provided 
here are thus unaffected by the actual shape of the temper- 
ature profile. 

3.2. Alternative scenario 
3.2.1. Balmer jump 
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Fig. 3. Plot of error bars for the pressure and the Balmer 
increment (AAZ^Baimcr)- This shows that a sharp increase 
in absorption depths at the limit of the Balmer series is not 
significantly detected. The increase in absorption depth of 
~0.030% is more likely to be due to Rayleigh scattering 
by gas at ^33 mbar pressure than by a Balmer jump of 
-0.030%. 

The rise of absorption depth at short wavelengths has al- 
ready been interpreted in terms of an extended cloud of ex- 
cited hydrogen that is optically thin at wavelengths shorter 
than the Balmer limit (Ballester et al. 2007). To analyze 
this possibility, we add one free parameter to the model 
with a jump in absorption depth between 3600 and 3900 A, 
Aj4Z)Baimer- The additional absorption is taken to be zero 
above 3900 A and linearly increases to reach A ADBaimor at 
3600 A and below. A fit to the data using this additional 
free parameter provides AAi:>Baimcr= 0.002%±0.007%. The 
best fit corresponds to a x'^ of 40, while a Balmer increment 
of '--^0.030% as claimed by Ballester et al. (2007) increases 
the to ^ 50 and decreases the pressure at zero altitude 
to ~15 millibars to avoid absorption by Rayleigh scattering 
(Fig. [3]). We therefore conclude that Rayleigh scattering 
with a smooth increase in absorption depth at short wave- 
lengths better fit the observational data than the "Balmer 
jump" with its sharp increase in absorption depth below 
- 4000 A. 

Alternatively, the "Balmer jump" in absorption due 
to a cloud of excited hydrogen could also be consistent 
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with the measurements, if we consider an additional ab- 
sorber to explain the low variation in absorption depth be- 
tween 4000 and 5000 A. Therefore, we add a second free 
parameter to the model, namely a minimum absorption 
depth, AD^in, which is expected when an optically thick 
cloud layer is present at high altitude. For these two ad- 
ditional parameters, we find that AL)min=1.4452±0.0016% 
and AAD^ 0.031%±0.007%. In this model, fit to the data 
again constrains the pressure at zero altitude to be lower 
than '^15 millibars. Since the sodium line profile is pres- 
sure broadened, a lower pressure at reference altitude re- 
quires a far higher sodium abundance to fit the blue tail 
of the absorption line profile. Using this lower limit for 
the pressure, we find that the sodium abundance must be 
X(Nal) = 1.5 X 10~^ which is 10 times the solar abun- 
dance. In addition, the lower pressure at zero altitude also 
constrains the pressure at which the core of the sodium 
line is detected to be below ~10^^bar; this is a very low 
pressure at which all atomic sodium should have disap- 
peared through ionization by UV flux, even at lower al- 
titude (Fortney et al. 2003). A fit to the data using a cloud 
layer and a "Balmer jump" therefore requires a high sodium 
abundance at significant height in the atmosphere where 
sodium is believed to be ionized. The "Balmer jump" sce- 
nario thus appears to be unlikely. 

We conclude that Rayleigh scattering is probably the 
simplest and most robust interpretation of the data because 
the absorption-depth curve from 3000 to 6000 A, including 
the rise between 4000 and 5000 A, can be interpreted by 
means of Rayleigh scattering and a sodium abundance close 
to solar. In contrast, the "Balmer jump" would require two 
additional components and a high sodium abundance. 

3.2.2. Atomic lines 
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Fig. 4. A similar plot to Fig. [l] with fit including absorp- 
tion by atomic iron lines. When the atomic iron lines are in- 
cluded, a satisfactory fit to the data can be obtained assum- 
ing a lower pressure of ~15 millibars and larger abundance 
for sodium and iron of about 10 times solar abundances. 
The dotted line shows the absorption depth obtained ne- 
glecting the Rayleigh scattering. This demonstrates that 
even with large iron abundance, the Rayleigh scattering 
has a significant contribution to the absorption depth level 
between 4000 and 5000 A, imposing a minimum pressure at 
the corresponding altitude. 



and magnesium, iron is supposed to be trapped in hydro- 
gen compound such as FeH. Because iron has the highest 
condensation temperature, it is not expected to be present 
in the gas phase in the absence of other species. 



We considered another alternative to Rayleigh scattering 
by the addition of plausible atomic lines to interpret the 
increase in absorption depth at the shortest wavelengths be- 
low 5000 A (Barman 2007) . Taking into account the species 
abundances, the strongest atomic lines in the 3000-5000 A 
wavelength range are those of Cai, Fei, All, and Cr i. Since 
the strong calcium fine at 4230 A is not observed, the abun- 
dance of atomic calcium must be extremely low. The same 
conclusion applies to aluminium at 3950 A and to chromium 
at 4270 and 3600 A. These species cannot be responsible for 
the observed absorption at short wavelengths. 

Nonetheless, atomic iron presents a series of strong lines 
from 2900 to 3860 A. When the atomic iron lines are in- 
cluded, a satisfactory fit to the data is achieved when we 
assume (i) a lower pressure of '^15 millibars; (ii) a higher 
abundance of both sodium and iron, of about 10 times so- 
lar abundance, and (iii) very low calcium, aluminium, and 
chromium abundances (Fig. [¥]). The higher sodium abun- 
dance is required to fit the blue wing of the line profile, 
because the line is pressure-broadened, and the pressure is 
lower. 

Calculating the best fit solution with and without iron 
lines (Fig. [1] and H]), we find similar differences between the 
data and fit with x^=40 for 52 degrees of freedom in both 
cases. This does not favor either one of the two models. 
However, the model without lines of atomic iron is pre- 
ferred because it does not require the large ten times solar 
abundance for sodium and iron. As for calcium, chromium 



4. Discussion 

It is noteworthy that the temperature is found to be rel- 
atively high compared to first estimates obtained from 
secondary transits measurements (Deming et al. 2005). 
However, a detailed model is needed to interpret the es- 
timated brightness temperature obtained from secondary 
transits and to determine the pressure and altitude at which 
this temperature corresponds (Burrows et al. 2005). The 
temperature obtained here, mainly constrained by the slope 
of the Rayleigh scattering, is consistent with temperature- 
pressure profiles obtained by some detailed models of the 
atmosphere of HD 209458b which provides high temper- 
atures up to 2000 K at 100 millibars, in the night side 
(Burrows et al. 2006). However, other models to interpret 
secondary transits measurements provide lower tempera- 
tures (Seager et al. 2005) and lower temperature gradients 
in the millibars pressure domain (Iro et al. 2005; Fortney 
et al. 2006). Nonetheless, a wide range of models are con- 
sistent with the data, and although efforts were made to 
consider the widest range possible (Seager et al. 2005), mea- 
surements remain limited and their interpretation is model- 
dependent, despite significant recent progress (Knutson et 
al. 2007; Burrows et al. 2007a). 

As shown above, the temperature estimates, in the 
Rayleigh regime, relies mainly on the slope and, there- 
fore, on the absorption depth measured in the 3000-4000 A 
range, where noise is the highest. Nevertheless, a fit to only 
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measurements for absorption depth above 4000 A provides 
also a high temperature of Tq = 2350 ± 300 K (1-cr) at 
zero altitude. Therefore, although the temperature given in 
Sect. [3] must be considered in hght of its large associated 
error bars, a high temperature is also required to explain 
the variation in absorption depth between 4000 and 5000 A 
where error bars are the lowest. In addition, this is con- 
sistent with the temperature required to explain the high 
sodium abundance below the sodium condensation level 
(Sing et al. 2008b). 

This fit, which uses only measurements above 4000 A, 
strengthens the Rayleigh interpretation. The variation in 
absorption depth between 4000 and 5000 A cannot be ex- 
plained by the absorption in the Balmer series, while the 
extrapolation of the 4000-5000 A measurements assuming 
Rayleigh scattering perfectly matches the absorption depth 
measured between 3000 and 4000 A. 

Importantly, even if Rayleigh scattering is not detected 
and the observed increase in absorption depth is due to 
another phenomenon, such as the Balmer jump (Ballester 
et al. 2007) or atomic lines (Barman 2007), the present 
calculations provide a strong upper limit to the pressure at 
the absorption depth level of 1.453%. The result of Sect. [3] 
shows that at this level, the pressure must be equal to or 
smaller than 33 millibars, a pressure at which because of 
Rayleigh scattering a grazing line of sight becomes optically 
thick at 4300 A. 

The red wing of the magnesium line at 2850 A could 
also explain the increase in absorption depth at the shortest 
wavelengths, however the abundances of both sodium and 
magnesium would have to exceed approximately '^20 times 
solar. We therefore conclude that the increase in absorp- 
tion depth at the shortest wavelengths is best explained by 
Rayleigh scattering. This represents the first direct detec- 
tion of the main constituent of the planetary atmosphere: 
molecular hydrogen. 

5. Conclusion 

Before, the determination of the planetary radius-pressure 
relationship was missing although this is needed to deter- 
mine absolute abundances. Until now, many models have 
relied on the assumption that pressure at the mean planet 
transit radius is close to 1 bar. Using Rayleigh scattering, we 
now determine that this pressure should be ~30 millibars. 
Even if the Rayleigh scattering does not dominate the ab- 
sorption, the pressure must be lower than ~30 millibars. 
Consequently, the radius of HD 209458b at 1 bar is deter- 
mined to be about 1.29 Jupiter radius, slightly lower than 
previously assumed. However, the difference is not sufficient 
to explain the puzzling large radius of HD 209458b (Guillot 
et al. 2006; Burrows et al. 2007b; Chabrier et al. 2007). 

In the near future, the new generation of instruments 
with increased capabilities in the near UV should enable the 
physical properties of a constantly increasing list of tran- 
siting extrasolar planets to be determined. When Rayleigh 
scattering by H2 in these planets is detected, it will provide 
the reference baseline for the determination of the absolute 
abundances of all the elements to be detected using absorp- 
tion spectroscopy (Ehrenreich et al. 2006). Measurement of 
the atmospheric composition will become an increasingly 
important tool to help understand the origin and nature of 
extrasolar planets, as increasingly smaller extrasolar plan- 
ets are detected (Gillon et al. 2007). 
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